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Abstract 

The localization kinetics of a regular block-copolymer of total length A'^ and block size M at a 
selective liquid-liquid interface is studied in the limit of strong segregation between hydrophobic 
and polar segments in the chain. We propose a simple analytic theory based on scaling arguments 
which describes the relaxation of the initial coil into a flat-shaped layer for the cases of both Rouse 
and Zimm dynamics. For Rouse dynamics the characteristic times for attaining equilibrium values 
of the gyration radius components perpendicular and parallel to the interface are predicted to scale 
with block length M and chain length A*' as oc M^'^'^'^ (here v ~ 0.6 is the Flory exponent) 
and as cx A'^^ , although initially the characteristic coil flattening time is predicted to scale with 
block size as oc M. Since typically N ^ M for multiblock copolymers, our results suggest that the 
flattening dynamics proceeds faster perpendicular rather than parallel to the interface, in contrast 
to the case of Zimm dynamics where the two components relax with comparable rate, and proceed 
considerably slower than in the Rouse case. 

We also demonstrate that, in the case of Rouse dynamics, these scaling predictions agree well 
with the results of Monte Carlo simulations of the localization dynamics. A comparison to the 
localization dynamics of random copolymers is also carried out. 
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1 Introduction 



The behavior of hydrophobic - polar (HP) copolymers at a selective penetrable interface (the interface 
which divides two immiscible liquids, like water and oil, each of them being favored by one of the two 
types of monomers) is of great importance in the chemical physics of polymers. For strongly selective 
interfaces, when the energy gain for a monomer in the favored solvent is large, the hydrophobic and 
polar blocks of a copolymer chain try to stay on different sides of the interface leading thus to a 
major reduction of the interfacial tension between the immiscible liquids or melts which has important 
technological applications, e.g. for compatilizers, thickeners or emulsifiers. Not surprisingly, during 
the last two decades the problem has gained a lot of attention from experiment [1-4], theory [5-9] as 
well as from computer experiment [10-14] . While in earlier studies attention has been mostly focused 
on diblock copolymers [2,4] due to their relatively simple structure, the scientific interest shifted later 
to random HP-copolymers at penetrable interfaces [6-9, 14]. Until recently though, the properties of 
regular multiblock copolymers, especially with emphasis on their dependence on block length M, have 
remained largely unexplored. This has been a driving force for our research, since experimentalists are 
now able to synthesize polymers with an increased variety of controlled structures. 

In a recent study [15], we showed that the equilibrium properties (structure, diffusion coefficient, 
etc.) of a regular HP-copolymer at a selective liquid-liquid interface are well described by a scaling 
theory in terms of the total copolymer length N (the number of repeating units in the chain) and 
the block size M (the number of consecutive monomers of the same kind) as well as the selectivity 
parameter x, that is, the energy gained by a monomer when being in the more favorable solvent. 

In terms of the selectivity parameter, there are three adsorption regimes that can be distinguished: 

• if X is smaller than a critical value Xc, the interface is too weak to affect the polymer at all so 
that the macromolecule conformation is identical to that in the absence of an interface; 

• for X > Xc the polymer starts to feel the presence of the interface so that the coil radius per- 
pendicular to the interface starts to shrink while parallel to the interface this radius swells. In 
this regime the scaling approach developed in [15] correctly describes the statics of the problem. 
This regime is labeled as of weak localization; 

• if X is increased furthermore, a second threshold value Xoo is reached. This value corresponds to 
an interface that is strong enough to induce a perfect flattening of the copolymer in which all the 
monomers are in their preferred environment. In this regime of strong localization, an increase in 
X does not change the shape of the polymer anymore since the polymer is already in equilibrium, 
as flat as it can get, depending on its structure, mainly on the length of the blocks. 
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Using scaling arguments, we demonstrated [15] that: (i) the crossover selectivity decreases with 
growing block length as Xc oc M~'^^^'^^^'^ (where v « 0.6 is the Flory exponent in three dimensions) and 
the threshold to the strong localization regime vanishes as Xoo M~^, and (ii) the size of the copolymer 
varies in the weak localization regime as i?g_L oc M"''(i+'')/(i"'') and Rg\\ oc M(''2"'')(i+'')/(1"'^), whereas 
Rg±_ OC , Rg\^ (X M''-"^'"'' in the case of strong locahzation (where 1/2 = 3/4 is the Flory exponent 
in two dimensions). We shall use these results in the present paper in which we suggest a theory of 
the localization kinetics of a regular block copolymer at a liquid-liquid interface based on dynamical 
scaling arguments. To the best of our knowledge, so far there have been no attempts to treat this 
problem analytically or by means of computer simulation. 

2 Scaling analysis 

In our dynamical scaling analysis wc consider a coarse-grained model of a multiblock copolymer con- 
sisting of N repeat units which is built up from a sequence of H- and P - blocks each of length 
M. For simplicity one may take the interface with negligible thickness as a flat plane which sepa- 
rates the two selective immiscible solvents. The energy gain of each repeat unit is thus — x, provided 
it stays in its preferred solvent, and the system is considered in the strong localization limit where 
X > Xoo oc [15, 19]. In our consideration, we skip the time needed by the polymer to move by 

diffusion through the bulk into the vicinity of the interface, and place the center of mass of an unper- 
turbed coil at time t — at the interface whereby the localization field is switched on. The initial coil 
will then start relaxing with time into a flat ("pancake") equilibrium configuration in-plane with the 
interface and the kinetics of relaxation will be determined by the sum of the various forces acting on 
the copolymer. 

It might be argued that a more ptiysical initial configuration could be represented by a polymer 
that, after diffusing in the bulk of one of the solvents, gets to the interface and touches it with one 
single monomer. We have indeed tried this initial configuration as well, and concluded that, once the 
adsorption process starts, it proceeds exhibiting the same physics (e.g. the same scaling behavior, see 
Section 4.1). The choice of our initial condition just eliminates the spread in the initial time at which 
adsorption sets on. This spread is caused by polymers which occasionally drift away from the interface, 
provided the coil touches the interface initially with the "wrong" kind of P- or H-monomers. 

2.1 Rouse dynamics 

In order to estimate the driving force of the flattening process, one may recall [15, 19] that the effective 
attractive energy e (per diblock, e.g. a pair of consecutive blocks of H- and P- monomers) in the strong 
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localization limit is e w xM- this case the diblock plays the role of a blob and the overall attractive 
free energy Fattr oc eM cx where N ~ N/M is the total number of blobs. Thus the effective 
driving force perpendicular to the interface is /^tr ~ ~XooN/ R± where R± denotes the perpendicular 
component of the radius of gyration. This force is opposed by a force of confinement due to the 
deformation of the self-avoiding chain into a layer of thickness R±. The corresponding free energy of 
deformation is simply estimated as i^conf — N{b/ R^Y^'^ where b is the Kuhn segment size [20]. For 
the respective force then one gets f^^i — Nb^^'^ / R^/'^^^ . 

The equation of motion for R± follows from the condition that the friction force which the chain 
experiences during the motion in the direction perpendicular to the interface is balanced by the sum 
of /ai;tr 9-^^ /conf • -'^^ ^hc casc of Rousc dynamics, each chain segment experiences independent Stokes 
friction so that the resulting equation of motion has the form 

where Co is the friction coefficient per segment. 

During the flattening process, the chain spreads parallel to the interface due to the excluded volume 
interaction. Within Flory mean - field arguments, the corresponding free energy F^^ — wiV^/(i?|i?j^), 
where v is the second virial coefficient and is the gyration radius component parallel to the interface. 
The corresponding driving force is /it — f iV^/(i?|i?j^). This term is counterbalanced by the elastic 
force of chain deformation, i.e. by /j^j ~ —R\\/{b'^N). 

Taking into account the balance of forces in the parallel direction, the equation of motion for R\\ 
takes then the following form 

Evidently, the excluded volume interactions provide a coupling between the relaxation perpendicular 
and parallel to the interface. Thus, eqs. Q and (0) describe the relaxation kinetics of a multiblock 
copolymer conformation at a selective liquid-liquid interface. One may readily verify that the equilib- 
rium solutions which follow from these equations are 

R°l ~ bM\ (3) 

and 

R"^^ :^{vbY'^M-^I^N^\ (4) 

These coincide with the equilibrium expressions for R± and R\\ derived earlier from purely scaling- 
based considerations [15]. The only difference is in the power of the M - dependence in eq. Q which 
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looks like M~^/^ instead of M"'^'^'^^^^ in ref. [15], albeit numerically the values of both exponents 
coincide: v/A ~ {v2 — i^) ~ 0.15. To get the full solution of the equations of motion it is convenient to 
rescale the variables as a; = R±/ {bM'^), y = R^^M'^ / (vbN^) so that eqs. and jSJ can be written 
in the dimensionless form 

dx 1 1 , 

The characteristic times for relaxation perpendicular and parallel to the interface in eqs. |(SJ - lO 
should then scale as rj_ ~ Co^^A^^^^", and t\\ ~ C^b'^N'^. 

Eq. is independent of y so that it can be solved exactly: 

^ [l-F{2v,l-l + 2v-x''''{t)) 

-^[l-^^(2^,l;l + 2^;xi/-(0))l =-- , (7) 

where F{a, /3; 7; z) is the hypergeometric function and x(0) — R±{0)/bM'^ is the initial value. In the 
early stages of relaxation (i.e. aX t <^t±) one has a; ^ 1 and the solution, eq. {Tj), with R±{0) ~ bN'^ 
reduces to 

Rl{t)c^Rl{0)--±- , (8) 
Co A/ 

so that the perpendicular collapse of the chain should last proportionally to the block length M. 

In the opposite limit of late stages kinetics, t t±^ and x > 1, that is, close to equilibrium, the 
relaxation of Rg{t) is essentially exponential with t±_ oc A/^-^: 

^-l + exp(-^) . (9) 

Moreover, as typically N ^ M, one expects that r|| ^ t±, i.e. the chain coil collapses first in the 
perpendicular direction to its equilibrium value and after that slowly extends in the parallel direction. 
In this case one can use x{t) ~ = 1 in eq.© and derive the resulting solution for the parallel 
component of Rg as 

y(t)-l = [y(0)-l]exp(-i/r||) , (10) 

where the initial value y{0) < 1. As one may see in Fig. ^ the time evolution of R± and i?|| which 
follows from the solution of the system of differential equations lO - resembles qualitatively the 
simulation data even though one should bear in mind that the time axis for the former is in arbitrary 
units (scaling always holds up to a prefactor) whereas in the simulation time is measured in MC steps 
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(MCS) per monomer (i.e., a MCS has passed after all monomers have been allowed to perform a move 
at random). The relaxation of looks somewhat slower and, as our data for longer chains show, this 
effect becomes much more pronounced with growing number of blocks. 

Concerning Fig. ^ note that there are several approximations on which our analytical derivations 
rely. First of all, our expressions for the attractive free energy, i^attr, and the confinement free energy, 
^conf , are only appropriate in the limit of strong flattening which occurs during the late stages of the 
localization process. Moreover, the free energy i^cv, due to excluded volume interactions, is somewhat 
underestimated if one uses the second virial term approximation only. That explains why both analytic 
curves in Fig. ^ qualitatively agree well with the simulation data at the late stages of the relaxation 
process whereas, not surprisingly, this agreement is missing during its initial stages. 



2.2 Zimm dynamics 

In the case of Zimm dynamics, the hydrodynamic long-ranged interactions couple the motion of the 
solvent to that of the coil within the framework of the non-draining coil model [16]. When the coil 
flattens, the solvent velocity u(r) follows the velocity of the coil segments (the solvent entrainment), 
as depicted schematically in Fig. [5] In order to calculate the friction forces which act during the coil 
flattening, we recall that the energy dissipation in an incompressible fluid is [17]: 
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rio d r{VkUi{r)+ViUkir)) 



(11) 



where u(r) is the solvent velocity field, r/o ^ Co/b is the viscosity, and the integral is taken over the 
interior of the coil. In our case the energy dissipation can be estimated as 
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where we have taken into account that typically the velocity gradients are V^Ua; = VzUj, ~ {dR\\ /dt) / R± 
and WxUz = 'S/yUz ~ {dR±/dt)/R\\, and the coil volume is The relevant friction forces follow 

from eq. 112|) by differentiation with respect to the corresponding velocities [18] 
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where C is a constant whose value will be found later from the initial conditions. Note that the ratio 
of the friction forces is proportional to the aspect ratio of the flattening coil. Since the friction forces 
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are now proportional to one another, in contrast to the case of Rouse dynamics, one of the equations 
of motion becomes redundant. Thus, after equating the friction and driving forces one derives a single 
relationship between the perpendicular and the parallel components of the gyration radius, i.e. 



Riit) = 



' R± 



+ 



c.l/"' + l 



MR, 



(15) 



C 



As usual it is more convenient to use the rescaled variables x and y, defined above. In terms of 
these variables, ea. ((TB|l becomes 

N'^C ( 1 



y 



1 



(16) 



The constant C in ca. l|15(l can be fixed by the initial conditions, i.e. y = yo provided that a; = xq, so 
that C ^xo (xo - 1/yo) M^+'^''/[(^l - l/xl^"^ N'^], and eq. ^ takes the form 

_i _ xo {xo ~ yo^) (1 - x-^^") 



(17) 



As a consequence, only one equation of motion is required for R±, while can then be obtained 
directly from eq. (|15|l . The equation of motion for R± can be derived, if one notes that the friction force 
H1H|I is balanced by the sum of the forces corresponding to the free energy contributions corresponding 
to the attraction fattr = ~XocN/ Rj_, and to the confinement, fconf = Nb^^'^ / R^/'^^^ . One may exclude 
dR\\/dt from eq. H13|l by taking into account the solvent incompressibility [17], namely V • u(r) = 0, 
with V^Wa; — ^yUy ~ {dR\\/ dt) / R\\, and VzU^ ~ ~{dR±/dt)/ R±. This yields a simple relation 
between the two velocities and the coil aspect ratio 
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Using this relationship in eq. ill^ . one obtains the equation of motion of R± as: 

dRi_ 1 h^l" 
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Making use again of the reduced variables x and y, eq. (|19|l can be recast into the form 
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(20) 



The resulting equation of motion l|2U|l . however, in contrast to the case of Rouse dynamics, provides 
no clear-cut scaling relationship for the typical relaxation time Tzimm with respect to N and M . Making 
use of the expression for y in eq. (|17|l and taking the starting conformation as that of an unperturbed 
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self-avoiding coil i.e. R±{0) = i?||(0) = bN" , [that is, xq = {N/MY and yo = b'^N'^'' M" / {vN^)], 
eq. (|20|l can easily be solved numerically for arbitrary chain length N and copolymer block size M . In 
Fig. ^results are shown for N ~ 128, M ~ 32, and compared to the predicted localization kinetics for 
the case of Rouse dynamics, eqs. (jSJ and 

Two major conclusions are suggested by a close inspection of Fig. |21 Evidently, the presence 
of hydrodynamic interactions in the case of Zimm dynamics tends to erase the asymmetry in the 
localization rates perpendicular and parallel to the interface. This is due to the strong coupling between 
perpendicular and horizontal degrees of freedom induced by the solvent. In addition, compared to the 
case of Rouse dynamics, the localization in the presence of hydrodynamic effects takes considerably 
longer, at least one order of magnitude, under otherwise identical conditions. The physical reason for 
this slowing down is the presence of high velocity gradients which develop in the incompressible solvent 
during the coil flattening. 

3 The Simulation Model 

For the Monte Carlo simulations we use an off-lattice bead-spring model that has been employed 
previously for simulations of polymers both in the bulk and near confining surfaces [21]. Recently, it 
was applied by us to the study of the static properties of block copolymers at a selective interface [15], 
therefore we describe here the salient features only. 

Each polymer chain contains N effective monomers connected by anharmonic springs described by 
the finitely extendible nonlinear elastic (FENE) potential, 

Ufene = -§R' In [1 - ^^^^] ■ (21) 

Here £ is the length of an effective bond such that Imin < i < £maxi with trnin — 0.4, Imax = 1 being 
the unit of length, and has the equilibrium value £o — 0-7, while R = imax ~ = ~ ^min — 0-3, 
and the spring constant K is taken as K/kBT = 40. The nonbonded interactions between the effective 
monomers are described by the Morse potential, 

Um = eM{exp[~2a(r - r„i„)] - 2 exp[-a(r - r™„)]} (22) 

where r is the distance between the beads, and the parameters are chosen as rmin = 0.8, ei\i — 1, and 
a — 24. Owing to the large value of the latter constant, U]\i{r) decays to zero very rapidly for r > r„iin, 
and is completely negligible for distances larger than unity. This choice of parameters is useful from 
a computational point of view, since it allows us to use a very efficient link-cell algorithm [21]. From 
a physical point of view, the interactions Ufene and C/m, eqs. (|21ll and (|22|l make sense when one 
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interprets the effective bond as a Kuhn segment, comprising a number of chemical monomers along the 
chain, and thus the length unit imax = 1 corresponds physically rather to 1 nm than to the length of a 



with the localization of a copolymer in good solvent conditions, in U m{t) we retain the repulsive branch 
of the Morse potential only by setting U]\i{r) = for r > Vmin and shifting UMi^) upwards by eM- 
The interface potential is taken simply as a step function with an amplitude x, 



where the interface plane is fixed at 2; = 0, and cr(n) = ±1 denotes a "spin" variable which distinguishes 
between P- and H- monomers. In studying the flattening dynamics of these chains, we always start 
with a configuration which has been equilibrated in a solvent that is good for both types of monomers. 
At time t — 0, the interface is switched on (so that it goes through the center of mass of the chain) 
and the flattening dynamics is then followed until the chain reaches its new equilibrium configuration 
at the interface. As explained above, this represents a good initial condition that makes the simulation 
of the flattening process quite efficient. Since we are interested in the behavior of these chains in the 
limit of strong localization, we have chosen x — where Xc is the crossover selectivity as obtained 
in Ref. [15]. We use periodic boundary conditions in the plane of the interface while there are rigid 
walls in the z - direction, where the simulation box extends from z = —32 to z = 32. Typically we 
studied chains with lengths 32 < < 512 and block lengths 1 < Af < A^/8 whereby all measurements 
have been averaged over 1024 different and independent equilibrated starting configurations. 

4 Simulation results 

4.1 The characteristic times of locahzation 

The results of the Monte Carlo simulations have been compared to the predictions of the Rouse model 
introduced above. It is not possible to compare the simulation results to the predictions of the Zimm 
model since these simulations do not take into account the hydrodynamic contributions that are so 
important in the Zimm treatment. 

Fig. 2] shows two typical snapshots for a chain with A^ = 128 and M = 8 on the verge of the 
localization threshold (x = 0.25, 77 = 3.48, where 77 = xM^'^+^^/^N^'^^"''^^ is the "natural" scaling 
variable of the problem [15]) and in the strong (x — 10,77 = 139.3) localization limit. The value of the 
crossover selectivity for this chain is Xc = 0.67 [15]. 



covalent C — C bond (which would only be about 1.5 A). Since in the present study we are concerned 




(23) 
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Figure El shows the typical dependence of R^{t) on time in the case of = 256, and various M 
values. These results have been used to determine the dependence of the flattening dynamics on the 
block length M. In order to determine the early stages dynamics, the initial slope of these curves has 
been obtained and compared with the theoretical prediction of eq. 0- The results are presented in 
Fig. El 

From Fig. EJi one may readily verify that the initial collapse of R±{t ^ t±) with time closely 
follows the predicted rate oc according to eq. © in the limit of N,AI 3> 1. Moreover, for 

sufficiently large block size M the A^-dependence of this initial rate disappears so that asymptotically 
this initial perpendicular collapse is governed by the length M only, as predicted by the theory. The 
disagreement between the prediction of the scaling theory and the results of the simulations for small 
values of polymer length and block length is not surprising since the predictions of the scaling theory 
are valid in the limit of — ^ oo and Af — > oo, with N/M large. So, we can expect the results of 
the simulations to have a better agreement with the predictions of the theory in the case of very long 
polymers with many long blocks. 

We can go one step further in the analysis of the dynamics of the perpendicular component of the 
radius of gyration and assume that the scaling law for the early stages collapse rate has the form 

. i M \ 

^^=^1^ TTT , (24) 



J4 , 

where rj^ oc N"^ and M* oc N''^ and the scaling function G{x) has the form 

{1 , x < 1 
(25) 
X , a: > 1 

According to our prediction, for M > M* the A^ dependence is dropped out, so that we come to the 
conclusion that ai =71. It is possible to obtain the A^ - dependence of t± at small M from Fig. Et-, 
and find ai — 0.57. When this is known, 71 is fixed and it is possible to replot Fig. EJi in terms of 
T^/rj^ and M/M* as a master curve - cf. Fig. El). Certainly, this master plot is not perfect as far 
as both A^ and M are still too small for the asymptotic relation, eq. H24|l . to be well satisfied but it 
appears that for growing A'^ the different curves group on a single one, revealing a clear cut crossover 
region at Af/Af* w 0.3. 

From the curves shown in Fig. El it is also possible to obtain the characteristic time for the 
perpendicular component of the radius of gyration corresponding to the late stages of adsorption. 
To this end the fits of t±_ are taken only after the initially unperturbed coil has sufficiently relaxed, 
R±{t) / Rj_{0) < 1/e. The two vertical lines shown in Fig. El for example, indicate the time interval 
used to fit the late stages behavior for the polymer with total length of A^ = 256 and block length 
M = 8. 
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As shown in Fig. E^, during the late stages of locaHzation at the interface one observes the expected, 
cf. eq. ©, characteristic scahng of t± cx M^+^'^ which progressively improves as the asymptotic limit 
is approached. As expected, the initial dependence of t± on N is seen to vanish as M and N become 
sufficiently large, in agreement with eq.®. 

In Fig. [7^ we also demonstrate (see triangle-down symbols) that the use of initial polymer con- 
figurations which barely touch the interface (rather than having their center-of-mass at t = at the 
interface) does not qualitatively alter the simulational results while it greatly increases the statistical 
disarray of the sampling. As mentioned in the beginning of Section 2, this motivated us to choose in 
all our computer experiments coils whose center-of-mass is initially placed at the liquid-liquid interface 
before the selectivity of both solvents is switched on. 

The scaling approach introduced above can also be used for the results concerning the late stages 
of the localization process (Fig. W^)- In this case the scaling law reads: 



Again the condition that at M > M* the N- dependence of t± is dropped out leads to the link between 
72 and 02; in this case: 72 = 02/(1 + 2i^) = 0:2/2.2. We can obtain 02 = 1.52 from Fig. [7^, then 
determine the value of 72, and finally plot a master curve now in coordinates tj_/t^ versus M/M*. 
Fig. [7|d shows the master curve as obtained following this procedure. Eventually, in Fig. |Slwe display 
the measured scaling of the late stage characteristic time for the spreading of the chain in the plane of 
the interface, T||, with block size M for different chain lengths 32 < < 512. Apart from some scatter 
of data for too small TV — N/M, one nicely recovers the relationship r|| cx N'^ whereby, once again, the 
initial M-dependence gradually diminishes as A/" 3> 1. 

4.2 Rouse mode analysis 

An interesting question about the relationship between localization rate and copolymer composition 
concerns the mechanism which drives the initially unperturbed coil into a localized state and makes 
it "feel" the presence of a selective interface. A useful insight in this respect may be derived by the 
detailed Rouse mode analysis of the mode relaxation times. 

The Rouse modes are defined as the Fourier cosine transforms of the position vectors, rn, of the 
repeat unit n of a polymer chain. For the model of discrete polymer chain considered here they can 




(26) 



where rj^ (x N"^ and M* (x N^'^ and the scaling function T{x) has the following limits 




(27) 
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be written in terms of collective coordinates Xp as [22]: 



1 ^ 

Xp(i) = — ^r„(i)cos 



N 

n=l 



(n - l/2)p7r' 
N 



,p-0,l,2,...,iV-l (28) 



The factor 1/2 is related to the nature of the chain that has free ends. The time-correlation function 
of the Rouse modes is given by 

%,{t) = (Xp(t)X,(0)) (29) 

where the angular brackets denote taking an average over different moments of time t' separated by 
an interval t. For phantom Gaussian chains with no excluded volume interactions one gets 

$^--(t) = (Xp(0)X,(0))exp(^-^^, p=l,2,...,iV-l (30) 
with {p,q ^ 0) coefficients (Xp(0)Xq(0))gaMss — Spqb'^ /[8N sin^ {pn /2N)] and relaxation times 

^gauss _ /g-j^-, 

P 12kBTsm^{pTT/2N) 

so that the typical decay times for the first modes scale as r^""^'* cx iV^. The mode corresponding to 
p — describes the self diffusion of the chain. For polymers with excluded volume interactions no such 
analytic results exist but one can show [23] that Tp cx N^+^''. 

It is interesting to examine how the presence of an interface between the two solvents affects the 
leading modes dynamics. In the absence of an interface we find that in equilibrium the first Rouse 
modes (1 < p < 10) of a homopolymer chain in a good solvent decay according to expectations, 
Tp CX (1/^)1+2" CX and (Xp(O)Xg(O)) cx (1/^)1+2". This is shown in Fig. Elwhere the resuhs 

pertaining to one of the components of the modes are shown for a polymer of length N = 128 and 
block length M — 16. Due to the isotropy of the system, the behavior of the other two components, 
corresponding to the y and z directions is identical. In part (a) of the figure one can verify that the 
amplitude of the Rouse modes decays with the mode number as predicted. The discrepancy that occurs 
for large values of p is expected and it is due to the inability of a coarse-grained model, as the one 
used in this study, to capture properly the dynamics of the system on too short wavelengths, since the 
largest p values correspond t! o motions that occur on very short distances. Fig. Elb) demonstrates the 
orthogonality of the modes since the only (Xp(O)Xi(O)) that is non-zero is the one corresponding to 
p ~ 1. The results shown in the figure refer to the x component of the modes, but identical results are 
obtained if the other two components are plotted. In contrast, for copolymers localized at the liquid- 
liquid interface one finds a clear asymmetry in their behavior, strongly dependent on the direction that 
is considered. 
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Consider first the modes components that are in plane with the interface, x and y. The same 
quantities, that were plotted in Fig. |51for a homopolymer in the bulk, are displayed in Fig. 1101 for a 
copolymer adsorbed at the interface. Again one can see that the decay of the amplitude of the Rouse 
modes with mode number p still obeys a power-law. However, the polymers are not anymore 3D 
but rather 2D self- avoiding random walks. Evidently, the relevant Flory exponent in two dimensions 
U2 = 3/4 describes well the dependence of the mode amplitudes on the mode numbers as demonstrated 
by the dashed line with slope oc —(1 -f 2v2) ~ —2.5 in Fig. EH It is very interesting to notice how the 
modes corresponding to large values of p follow the power-law decay observed in the three-dimensional 
case. The modes with p < 8 describe monomer motions which involve more than M particles. Thus 
the modes corresponding to long wavelengths which span more than one block length feel the presence 
of the interface in contrast to those corresponding to shorter wavelengths (and show a dependence like 
cx see Fig. In part (b) of the same figure, we show that the Rouse modes components 

parallel to the interface are still orthogonal to each other. The only non-zero (Xp(O)Xi(O)) is the one 
corresponding to p ~ 1. This holds for both the x and y components of the modes. Careful inspection 
of the figure reveals that the value of {Xi{Q)Xi{Q)) is almost twice larger than that of the bulk case 
shown in Fig. |51 This is a consequence of the increased in-plane fluctuations that the polymer chain 
exhibits when it is confined to a plane. 

If we look at the component of the Rouse modes corresponding to the direction perpendicular to 
the interface, the results are quite different as shown in Fig. In Fig. Illf a^ one finds no power-law 
dependence of (Zp(O)Zp(O)) on the mode number p. The amplitudes of the perpendicular components 
of the modes decrease weakly with p and are several orders of magnitudes smaller than those of the 
parallel components of the modes. Evidently, the strong interface supresses all fluctuations of the 
chain perpendicular to the interface which is manifested by the strong decrease in the amplitude of 
the perpendicular components of the modes. 

Another striking difference in comparison to the results observed for the parallel components, 
can be seen in Fig. Illb . The modes are not orthogonal to each other anymore since the quantities 
{Zp{Q)Zi{Q)) , corresponding top — 3, 5, 7 up to p = 15, are non-zero and of the same order of magnitude 
as (Zi(O)Zi(O)). Thus one can conclude that the odd- numbered modes are coupled to one another. 
The amplitude of these couplings is so small, though, that the system does not lose its ergodicity. 
The same can be observed for the even-numbered modes too but the amplitude of these terms is even 
smaller than the amplitude observed for the odd-numbered modes. 

Another observation concerns the average values of the perpendicular components of the Rouse 
modes. If one calculates {Xp) or (1^), then due to the isotropy of the system in the plane of the 
interface these values vanish for all p's because there is no preferential orientation in the plane. If 
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one calculates {Zp), instead, the results are not trivial. The values obtained for a polymer of length 
N = 128 and block length M — 16 are shown in Fig. 112b . The squares represent the results obtained 
from the Monte-Carlo simulations of such chains. A pattern is clearly present in the (Zp): all the 
values corresponding to even p's vanish while those corresponding to odd p's are non zero. The values 
of (Zp) for odd p increase with growing p, and they change sign for p > N/M . In order to test if we 
were able to describe this behavior simply by the localization pattern at the interface, we approximated 
the shape of the localized chain perpendicular to the interface as a simple sine wave, 

■(n- l/2)7r" 



Z{n) = A{M) sin 



M 



(32) 



where n is an index that numbers the monomers in the chain [n = 1..N), and A is a prefactor related 
to the size of the adsorbed polymer in the direction perpendicular to the interface. This simplest 
sinusoidal shape, see Fig. 112b . is periodic with the block length AI while the first and last monomers 
of the chain are not pinned at z = accounting for the loose first and last blocks with respect to 
the interface. Using eq. (|32|l in eq. 128|) to obtain (Zp), we get what is represented by the circles 
in Fig. 112b . Note that we have rescaled the Monte Carlo results so that the two sets of data start 
from the same value for p = 1. Evidently, our simple expression for the average shape of the chain is 
able to reproduce remarkably well the simulation results, and in particular, to catch the trend in the 
amplitude growth as well as the change in sign observed at p = N/M. One may therefore conclude 
that the observed behavior of the Rouse modes for the chains adsorbed at the interface is indeed a 
consequence of the interplay between the interface and the regular block structure of the polymers 
under investigation. 

Eventually, we end our discussion of the Rouse-mode behavior with an analysis of the Rouse time 
spectrum - cf. Fig. ^| In both parts of Fig. E| we plot the relaxation times of the first ten Rouse 
modes for a polymer chain in the bulk (black circles in the figure). The dashed line represents the 
expected theoretical dependence of Tp on p as ~ p~^-^. As it can be easily observed, the results of the 
simulations follow the theoretical prediction quite nicely. As in the case of the amplitude of the Rouse 
modes, one observes deviations from the predicted law when the condition p/N <C 1 does not hold 
anymore. We compare the results for the chain in the bulk to the respective results for chains localized 
at the interface. In Fig. llSf a'). we consider the components parallel to the interface. Apparently the 
presence of the interface does not affect the relaxation time Tp considerably. Evidently, the times Tp 
for the chain at the interface are shorter than those for the chain in the bulk so that the presence of 
the interface effectively speeds up the inplane dynamics. In Fig. llSf b'). instead, the bulk behavior is 
compared to the behavior of the perpendicular component at the interface. In this case the times Tp 
corresponding to small values of p change more significantly. The presence of the interface accelerates 
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the relaxation of the modes that are characterized by large wavelengths (small values of p) since then 
the chain can better "sense" the interface. 

4.3 Localization of random copolymers 

Having understood the impact of block structure on copolymer localization kinetics, it is tempting to 
contrast it to the case of random HP-copolymers and the rate of their flattening onto a liquid-liquid 
interface. In recent years there has been a number of studies of random copolymers both by analytical 
theory [5,7-9] as well as by means of computer simulation [14,24,25]. As a rule, all these studies focus 
on the static conformational properties of random copolymers at the penetrable interfaces between 
two solvents. Thus it is now well established [5,24] that the perpendicular component of the gyration 
radius R±^ (x) of a random copolymer should scale with the degree of interface selectivity as a function 
of ?7 = for small and moderately large values of the interface selectivity x- 

In order to determine the crossover selectivity Xcr in our model, we carried out simulations on 
random copolymers with length 32 < < 512 and the results may be seen in the scaling plot of Fig. 
1141 where all data are found to collapse onto a single curve. One may estimate from the master curve 
in Fig. El that XcrV^ w 20. The dashed line in Fig. 1141 shows the predicted [5,24] dependence of 
the perpendicular component of the radius of gyration on the scaling variable cx rj^^" in the weak 
localization regime. As extensively explained elsewhere [15], the scaling is expected to fail when x is 
very large because then the blob size reaches its minimum possible value corresponding to a single 
block. Mathematically, the scaling fails for large values of the selectivity because of the different 
normalization factors used for the different curves as R±{Q) depends of course on A^. We have also 
checked another scaling hypothesis put forward for the first time in [24]: the total polymer density in 
the region of the interface should scale with when x li^s within the limits of the weak localization 
regime: 

pp{z) + Ph{z) = Pq9[{z - z,nt)x^''] (33) 

where pp{z) is the density of P monomers at a distance z — Zint from the interface, where Zint is the 
position of the interface {zint = in our simulations), ph{z) is the corresponding quantity for the 
H monomers, po is the total monomer density at the interface, and g{y) is a function close to being 
exponential. Fig. 1151 demonstrates that this scaling law is indeed observed for a wide range of polymer 
lengths and interface strength values. 

Also in the case of random copolymers, we performed all our studies of the localization kinetics 
in the limit of strong localization. It can be easily proved that, in the limit of infinite chain length, 
the average block size in a random copolymer is {M) ~ 2. In the case of polymers with finite length. 
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there are finite-size corrections that have to be taken into account [26]: the average block size is 
then always smaller than 2, albeit very close to this value, so it is appropriate to compare the rate 
of flattening of random copolymers to that of a multiblock copolymer with block size Af = 2 at the 
same selectivity strength x — '^Xcr- The results for the flattening dynamics of random copolymers are 
compared to those for regular block copolymers in FiguresElto|Hlwhere the results corresponding to the 
random copolymers are represented by filled symbols. In these figures it is possible to notice a marked 
difference between the early stage and the late stage dynamics. As seen in Fig. Et, the values of 
for the random copolymers and the corresponding block copolymers are very close for small values of 
N. As N increases, however, the difference between these two values grows as well and is particularly 
large in the limiting case of = 512 when the characteristic time for the random copolymer is about 
twice as large as the value for the corresponding block copolymer. A careful inspection of Fig. 
shows that the value of t± for a random copolymer with A^ = 512 is similar to that of a regular block 
copolymer of the same total length and a considconsiderablyr block size M — 8. 

The trend shown by the simulation results as the length of the polymers increases can be easily 
understood through a careful analysis of the nature of the random sequences that can be created in 
a polymer of finite length [26]: when the polymer length is short, e.g. A^ = 32, the probability of 
generating sequences of consecutive like monomers much longer than 2 is extremely small. As a matter 
of fact, the polymer will almost always contain a sequence of H— and P— monomers in which the 
length of each segment is small since every single polymer used to run the simulation is prepared in 
a neutral state with as many P— as i?— monomers. On the contrary, when a very long polymer is 
generated, the probability of having longer sequences of consecutive H— or P— monomers is increased. 
A long random copolymer will likely have long sequences of like monomers at some point along the 
chain. These long sequences are responsible for the much faster adsorption kinetics as compared to 
that of the regular block copolymers. The presence of long sequences of like monomers along the 
chain apparently affects only the early stages of the localization dynamics: from Figures □i andEJi 
it can be observed that during the late stages of the adsorption process, the behavior of random 
copolymers is almost identical to that of the corresponding regular block copolymers with M = 2. 
This observtiobservationunderstood by considering that during the late stages of adsorption the long 
blocks that are present in the random copolymers have already been adsorbed at the interface and 
it is the remaining short blocks that drive the system in order to complete the adsorption process. 
During this stage it is then the shortest blocks that are important for the adsorption process, and in 
this respect the random copolymers are very similar to the regular blocky ones, with blocks of length 
equal to 2. 
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5 Conclusions 



Prom the results presented above, it appears that the simple scaling theory that we have developed to 
describe the kinetics of copolymer adsorption on a selective liquid-liquid interface is able to capture 
quite well the most salient features of the problem as a comparison with extensive MC simulation data 
demonstrates. The agreement between the predictions of the model and the results of the simulations 
in both limits of early and late stages of adsorption is remarkable especially for long polymer and block 
lengths. 

The analysis that we have carried out regarding the Rouse modes of the system indicates that 
polymer chains that are strongly absorbed at the interface behave as two dimensional random walks 
in the plane of the interface, and show almost no fluctuations in the direction perpendicular to the 
interface. A simple trial function for the shape of the chain in the direction perpendicular to the 
interface is sufficient to capture the features of the modes that are specific to the adsorbed chains. 

Finally, we have also shown how the random copolymers behave at the interface. In particular, 
we have been able to prove how the presence of long blocks in the random sequences of the polymer 
chains affecs maff'ectsthe early stages of adsorption while the behavior of the random chains in the late 
stages of adsorption is almost identical to that of regular block copolymers with M = 2, as expected. 

Although the present study deals with the simplest situation of symmetric interface selectivity and 
copolymer composition and it assumes equal viscosity of both immiscible liquids, it can be viewed as 
a first step in exploring the localization kinetics of copolymers at the liquid-liquid interface. We also 
believe that our results could stimulate exciting laboratory experiments and be tested, for example, 
by means of the method of videomicroscopy of single macromolecules [27] . 
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Figure 1: Relaxation of R± and i?|| with time, following from Monte Carlo simulations of a copolymer 
of length N = 128 and block size M ~ 16. Both components are normalized by their respective 
equilibrium values. The inset shows the result for the same quantities (and the same set of parameters) 
as obtained from the numerical solution of eqs. (0) and jnj. 
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Figure 2: The solvent velocity profile in the Zimm model due to the cntrainmcnt of the solvent by 
the coil segments. The corresponding velocity gradients are V zU-x — ^ zUy — {dR\\/ dt) / R±^ and 
V^w^ = \JyU^ ~ {dR±/dt)/R\\. 
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Figure 3: Log- log plot of the normalized perpendicular and parallel components of the gyration radius 
vs time elapsed after the onset of localization for a regular block copolymer of length N = 128 and 
block size M = 16 at a flat liquid- liquid interface for the cases of Rouse and Zimm dynamics. Here 
V = and time is dimensionless, measured in units of (b'^. 
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Figure 4: Snapshots of typical configurations of a copolymer with N = 128, i\/ = 8 at x = 0.25 (a), 
and X = 10 (b). The value of the critical selectivity for this chain is Xc — 0.67. Here the light beads 
prefer the solvent that occupies the upper half of the simulation box, while the darker beads prefer the 
lower half of the box. 
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Figure 5: Typical dependence of R±{t) on time in the case of = 256, and M — 2 (solid line), 4 
(dotted line), 8 (dashed line) and 16 (dot-dash line). Arrows denote the time interval where t± for the 
late stage analysis has been determined by regression in the case of M = 8, see text for details. 
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Figure 6: (a) Variation of the characteristic time for the initial relaxation of R± with block length M 
for chains of length 32 < < 512. The prediction of the scaling theory is shown by a dashed line. 
The filled symbols represent results for random copolymers with the same values of N. (b) Master 
scaling plot for the localization of chains of different length N and block size M. 




Figure 7: (a) Variation of t± with block length AI for chains of length 32 < < 512. The slope 
of the dashed line is « 2.2, according to eq. 1^. The filled symbols represent results pertaining 
to random copolymers with the same chain length N. The triangle-down symbols represent results 
corresponding to initial configurations in which the chain barely touches the interface at t = 0. When 
not explicitly shown, the error bars are smaller than the size of the symbols, (b) Master scaling plot 
for the localization of chains of different length N and block size M. 
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Figure 8: Variation of T|| with chain length N for blocks of size M — 2, 4, 8 and 16. The slope of 
the dashed line is 2, as predicted by the theory. The filled triangles represent results for random 
copolymers with the same values of N. 



10 



<x > 

p 



10 



10 



(a) 



10 




Figure 9: Rouse modes for a chain in the bulk: (a) (Xp(O)Xp(O)) and (b) (Xp(O)Xi(O)) as a function 
of the mode number p for a chain with iV = 128 and M = 16 in the bulk of a solvent that is good for 
both kinds of monomers. The dashed line in (a) shows the predicted scaling law oc p~^'^. 



25 



10 



<x' > 10'' 
p 



10"' 



10"- 



(a) 



-2.5 



- -2.2 



10 



6 
5 

4- 



<X,X >3- 

1 p ^ 



2- 
1- 

Oo 1 



(b) 



Figure 10: Rouse modes for a chain at the interface: (a) {Xp{0)Xp{0)) and (b) (Xp(O)Xi(O)) as a 
function of the mode number p for the components of the Rouse modes in the direction parallel to the 
interface for a chain with A'' = 128 and M = 16 strongly adsorbed at a selective interface. The dashed 
lines in (a) show the predicted scaling laws. 
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Figure 11: Rouse modes for a chain at the interface: (a) {Zp{Q)Zp{0)) and (b) {Zp{0)Zi{0)) as a 
function of the mode number p for the components of the Rouse modes in the direction perpendicular 
to the interface for a chain with N = 128 and M = 16 strongly adsorbed at a selective interface. 
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Figure 12: (a) Comparison between the average values of the perpendicular components of the Rouse 
modes obtained from the trial function, see eq. H32() . and the values obtained from Monte Carlo 
simulations for a chain with N ~ 128 and M = 16 strongly adsorbed at a selective interface; (b) the 
trial function used to describe the average conformation of a chain with N — 128 and M = 16. 
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Figure 13: Relaxation times of the Rouse modes: (a) comparison between the bulk behaviour and the 
behaviour of the components parallel to the interface for the localized chains; (b) comparison between 
the bulk behaviour and the behaviour of the component perpendicular to the interface for the localized 
chains. As above, the results shown were obtained using a chain with N — 128 and AI — 16 strongly 
adsorbed at a selective interface. 
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Figure 14: Rj_{r]) / R±{0) as a function of the scaling variable 77 — xv^ for different chain lengths. 
The dashed line represents the predicted scaling behavior in the weak localization regime, oc rj~^'^ . 




Figure 15: Scaling plot of the total monomer density in terms of [z — Zint)~x^'^ for different combinations 
of polymer length and interface strength. 
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